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SUMMARY 

An Important problem that has emerged from combined analytical/experimental 
investigations is the task of identifying and quantifying the differences 
between results predicted by F.E. analysis and results obtained from experiment. 
The objective of this study is to extend and evaluate the procedure developed 
by Sidhu for correlation of linear F.E. and modal test data to include struc- 
tures with viscous damping. The desirability of developing this procedure is 
that the differences are identified in terms of physical mass, damping, and 
stiffness parameters instead of in terms of frequencies and modes shapes. 

Since the differences are computed in terms of physical parameters, locations 
of modeling problems can be directly identified in the F.E. model. 

From simulated data it was determined that the accuracy of the computed 
differences increases as the number of experimentally measured modes included 
in the calculations is increased. When the number of experimental modes is at 
least equal to the number of translational degrees of freedom in the F.E.- model 
both the location and magnitude of the differences can be computed very accu- 
rately. When the number of modes is less than this amount the location of the 
differences may be determined even though their magnitudes will be under 
estimated. 


INTRODUCTION 

The dynamic characteristics of structural systems are often predicted 
using Finite Element (F.E.) analysis and then later verified experimentally 
with dynamic analysis testing systems. Increased demands for reliability, 
minimal vibrations, optimum performance, and low cost design, among other cri- 
teria, have increased designers needs for sophisticated dynamic analysis test- 
ing techniques. Since the 1960‘s F.E. computer programs have become the 
preference of designers for analytical dynamic analysis. The use of F.E. com- 
puter codes has become especially widespread in the automotive and aerospace 
industries due to the requirement to analyze very large and complex structures. 
Commercial F.E. computer programs such as NASTRAN, ANSYS, and SAP (ref. 1) are 
available to anyone having access to a computer terminal. 

In many situations experimental verification is required to insure the 
validity of the results predicted by the F.E. analysis. Aerospace structures, 
which are very expensive and have rigorous safety and reliability requirements 
normally require experimental verification (ref. 2). Automobile prototypes 
are also experimentally verified to insure that vibration and noise problems 
will not exist in production models. Hundreds of other applications of F.E. 
analysis and experimental validation can be found. 






Digital signal analyzers are the most commonly used systems for experi- 
mental verification. Digital signal analyzers, which utilize the Fast Fourier 
Transform (FFT) developed In the 1960's (ref. 3), allow rapid and relatively 
accurate determination of structural transfer functions, resonant frequencies, 
and characteristic mode shapes. Modern digital analysis equipment has both 
automated the modal extraction process and decreased the required data acquisi- 
tion and post-processing time. These systems have replaced traditional analog 
devices because of their high speed and their ability to measure many modes 
simultaneously. 

An Important problem that has emerged from these combined analytical/ 
experimental Investigations Is the task of Identifying and quantifying the dif- 
ferences between results predicted by F.E. analysis and results obtained from 
the experiment. Although both the F.E. and experimental methods can be 
accurate from a theoretical standpoint. Inaccuracies do exist In their applica- 
tions to real structural problems. In the case of F.E. modeling there Is con- 
siderable uncertainty In the modeling of Items such as boundary conditions, 
joint flexibilities, and damping. Because of this, the F.E. results are not 
exact since the Input data Itself Is approximated. Also, It Is not possible to 
completely eliminate experimental error. F.E. analysts take the responsibility 
for producing theoretically correct computer codes but sometimes do not place 
enough emphasis on predicting the behavior of real world structures. The 
experimentalist, through testing, often show the limitations of the F.E. 
analysis, but do not always present clear cut procedures for quantifying the 
differences In a useful manner. 

A communication gap can exist between the experimentalist and the F.E. 
analyst when the experimentalist can not provide the quantitative data required 
by the analyst to Identify the differences between the experimental data and 
the F.E. model. The gap exists because the experimentalist normally measures 
frequencies and mode shapes In a vibration test, while the analyst requires a 
mass, damping, and stiffness matrix for describing the F.E. model. 

It would be useful If the differences between the experimental data and 
F.E. model could be found In terms of discrete mass, stiffness, and damping. 

If this could be done, and the experimental data was reliable, a more accurate 
F.E. model with Improved mass, damping, and stiffness descriptions could be 
created. This model could then be used for not only subsequent dynamic anal- 
ysis, but also for static analysis, for studying the effects of structural 
modifications, or for any analysis requiring the use of a mass, damping, or 
stiffness matrix. It would be Ideal If the discrete parameters could be 
measured experimentally but this Is not practical. For example, to measure 
the values for a row or column In the stiffness matrix, a displacement would 
have to be applied to the real structure while every other degree of freedom 
was constrained, and then the forces at all the other degrees of freedom would 
need to be measured. This would be both time consuming and require elaborate 
fixtures and Instrumentation. Experimental measurement of the mass and damping 
matrix would be at least equally complex. If not Impossible. 

One possible way to compare the experimental results to the F.E. model Is 
to compute analytical frequencies and mode shapes from the F.E. equation of 
motion and then compare them to the frequencies and mode shapes obtained from 
the experiment. The limitation of making a comparison at this level Is that 
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even though disagreements can be identified, the cause of the disagreements 
namely differences In the mass, damping, and stiffness matrices, can not be 
Identified of quantified. 

A more useful comparison between F.E. and experiment can be made through 
the equations of motion. By using the original F.E. equations and the equa- 
tions of motion derived from the experimental data, differences between experi- 
ment and F.E. coefficients can be Identified and corrected. Unfortunately, the 
procedure of deriving an equation of motion from the experimental frequencies 
and mode shapes Is not straight forward. To derive the equation of motion from 
experimental data requires that the same number of modes as degrees of freedom 
In the F.E. model be experimentally measured and that the experimental data not 
contain any measurement error or noise. If both of these requirements are not 
met the experimental data can not be used to construct a correct equation of 
motion. Since the coefficients for the equations of motion are computed by 
Inverting matrices containing the experimental mode shapes, these matrices must 
be square. In a typical experiment, the number of measured modes will not be 
equal to the number of degrees of freedom so the modal matrices will be 
rectangular Instead of square. Another difficulty Is that the experimental 
data will always contain some amount of experimental error and noise which 
makes the outcome of a matrix Inversion questionable. Also, If the highest 
modes In the structure are not Included In the experimental data the stiffness 
matrix computed from a modal matrix Inversion will be Incorrect (ref. 4). 
Finally, It Is difficult to measure the values of the mode shapes corresponding 
to every degree of freedom used In the F.E. model. This causes the order of 
the experimental matrices to be less than those In the F.E. equations. 

Previous research In this area has focused on using experimental data to 
Improve F.E. models rather than on Identifying the differences. Most of the 
techniques have been based on some form of a least squares fit. In the work by 
Berman and Flannelly (ref. 4), the analytical matrices are assumed to be close 
to the actual solution and then the smallest change In the analytical model 
that makes the experimental and analytical frequencies and mode shapes Iden- 
tical Is found. This assumption will not necessarily lead to an analytical 
model that Is physically representative of the actual structure. The only 
assurance is that the revised model will correctly predict the modes that were 
measured. The problems arising from using "Incomplete” data (data containing 
fewer modes than there are In the F.E. model) are also discussed In this work. 
In reference 5, Fuh, Chen, and Berman use similar approaches for correcting 
structures with viscous damping. 

Chen, Peretti, and Garba (ref. 6) refined a F.E. model of the Galileo 
spacecraft by first performing static tests to Improve the stiffness matrix, 
and then dynamic tests for correcting the mass matrix. The mass matrix cor- 
rection was based on a minimum change criteria. The limitations of this 
approach are that two Independent sets of tests must be run, and again, there 
Is no guarantee that actual physical characteristics will result from the least 
squares approach. 

Hart and Yao (ref. 7) discuss the advantages of using weighted least 
squares and Bayesian estimation. By using these extended forms of least 
squares methods, uncertainties In both the experimental data and analytical 
model can be Included In the updating procedure. It can be very Important to 
define the uncertainty in the experimental data since this data often contains 
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more error than the F.E. description. It does not make much sense to attempt 
to Improve a F.E. model with experimental data that Is less certain than the 
F.E. model. By Including relative uncertainties In the procedure, changes to 
the analytical model will not be applied Indiscriminately and the possibility, 
exists for retaining the physical meaning of the structure In the updated 
model . 

Dobb, Blakely, and Gundey (ref. 8), and Blakely and Walton (ref. 9) 
applied the Bayesian estimation procedure to a F.E. model of an offshore plat- 
form and a dam. In their study the effects of change In the uncertainties In 
both the experimental data and the F.E. parameters were Investigated. Unfor- 
tunately, well defined procedures do not exist for quantifying uncertainties 
so they had to be estimated using engineering Judgement. 

Sidhu (ref. 10) developed a procedure for approximating the difference 
between experimentally measured frequencies and mode shapes and F.E. parameters 
In terms of differences In mass, damping, and stiffness matrices. This 
approach has the potential for providing a direction to correct a F.E. model 
while retaining the physical characteristics of the real structure. The objec- 
tive of this work Is to extend the procedure developed by Sidhu for correlation 
of linear finite element and modal test data to Include structures with viscous 
damping. In this study, the derivation of the extended procedure and several 
case studies which use simulated experimental data are presented. The purpose 
of developing this procedure Is to formalize a process for Identifying the dif- 
ferences between experimentally measured frequencies and mode shapes and F.E. 
models In terms of differences In mass, damping, and stiffness. 


FORMULATION OF EQUATIONS 

The free vibration equation of motion for a damped, linear system can be 
written as: 


[M]{u} + [C]{u} + [K]{u} = {0} (1) 

where [M] Is the mass matrix, [C] Is the viscous damping matrix, [K] Is the 
stiffness matrix, and {u}, {u}, and {u} are the acceleration, velocity, and 
displacement vectors, respectively. The size of [M], [C], and [K] are nxn 
and {u}, {u}, and {u} are of size n, where n Is the number of degrees 
of freedom In the equations of motion. 


In only special cases can equation (1) be decoupled using normal modes 
(ref. 11). In general, when damping Is present, the solution of this equation 
results In complex eigenvalues and eigenvectors appearing In conjugate pairs. 
Since there are pairs of roots there will be twice as many roots as there are 
displacement degrees of freedom and the modal matrix will be of the order nx2n 
Instead of nxn. This rectangular modal matrix can not be used to decouple 
equation (1). Equation (1) can be rewritten In state vector form as: 


where {y} 


[A]{y} + [B]{y) = {0} 


({u}) [0] [M] 

\w) = [M] [C] 


[B] 


'-[M] [0]‘ 

[0] [K] 


( 2 ) 
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([A] and [B] are of order 2n x 2n and {y} 1s of order 2n.) 

The advantage of writing the equation of motion In state vector form Is 
that the modal matrix can now be used to decouple the equation. Assuming a 
solution {y} = {$}est and substituting into equation (2) leads to the 
eigenvalue problem: 


{[A]s + = {0} (3) 

For less than critical damping, the solution of this equation yields 2n com- 
plex eigenvalues s^, where Sp = -Wp Cr ^WQp. Wp Is the natural fre- 
quency, Wqp Is the damped natural frequency, and Cp Is the damping ratio for 
mode r. An equal number of complex eigenvectors are also obtained. 

Substituting the modal matrix [§] Into equation (3) and premultiplying 
by [§]T leads to: 


[§]f[A][§][\S'^l + [§]T[B][$] = {0} 


(4) 


from orthoginallty 

[$]T[A][$] = [^a^] and [$]T[B][$] = [\b>>^] 
where [^a->..] and ['^b>^] are diagonal matrices. 

If the eigenvectors are normalized with respect to the [A] matrix then: 

[§]T[A][$] = [I] • (5) 

and 

[§]T[B][$] = -[^s^] (6) 

Since the objective Is to determine the differences between the experi- 
mental model and the analytical model we need to find a common ground that 
will allow the comparison of the structural matrices computed from the F.E. 
analysis to the experimental frequencies and mode shapes. The differences 
between the F.E. ''[B]“ matrix and the [B] matrix computed from the experi- 
mental data (assuming that a [B] matrix can be created from the experimental 
data) Is written as: 


[D]b = [Blexp - [BIf.E. 


(7) 


rearranging 


[Blexp = [BIf.E. + [B]b 


then Inverting both sides 
and factoring out [B]p,E, 
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and 


[B] 


-1 

exp 


{[B]p {[I] + [B]p^p^ [D]g}} 


[Bi;Jp = an ^ <=> 

If the bracketed expression in equation (8) Is expanded using a Taylor 
series (refs. 12 and 13) and terms past the first derivative are dropped, 
equation (8) can be approximated by: 

[B]-' » ([II - 


multiplying out 


[B] 


-1 

F.E. 


- [B1^'.e,[B1b[B]-'.p, 

and then rearranging 

-tBlexp ^ » l 81 F^.t“lBtelF'.E. 


and solving for [DJg 

[DIb MB]p p ([B)^’e. - t<p't6lF.E. <1“> 

Using equation (6) to obtain [8]“^ and substituting Into equation (10) 
the final expression for the difference matrix [D]g Is obtained: 

[OIb » (B]p p a*lp p [^s^lp^J [*lj p - [»l„pt'-s-^i;'xp t»lIxpl[B'E.E. 

(11) 

The same approach can be used for deriving the difference In the [A] 
matrix. In this case: 

[01a == tAlF.E.<[‘lF.E.t‘lf.E. ' W*'exp't'')F.E. 

The format of equations (11) and (12) are well suited for computing the 
differences between the F.E. model and experimental data. Since these equa- 
tions do not require any Inversion of the modal matrices, the fact that all the 
modes are not measured does not cause a problem. As discussed previously, the 
modal matrix will not be completely known since fewer modes than degrees of 
freedom are typically measured. An Inversion of the frequency matrices are 
required, but this does not present any problems since these matrices are 
diagonal and their Inverses are just the reciprocal of the diagonal terms. 
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Once [0]f^ and [D]b are computed, the disagreement between the F.E. 
and experimental descriptions of the structure can be found. Since there Is a 
direct relationship between the elements of the mass, damping, and stiffness 
matrices and the elements of the [A] and [B] matrices, the discrepancies In 
mass, damping, and stiffness at any degree of freedom In the structure can be 
found by merely picking out values from the [D]/^ and [D]b matrices. For 
example, the disagreement In damping at the first degree of freedom would be 
obtained from the [D]/\ matrix at location [D(n + 1, n + 1)]a» the mass dis- 
agreement at [0(1, n + 1)1a. and the stiffness disagreement at [0(n + 1, 
n + 1)]b- Note that the mass discrepancy can be found from either one of 
two partitions In the [0]^ matrix or the [D]b matrix. 

It was mentioned previously that In preactive experimental mode shape data 
will normally not be available at all of the degrees of freedom used In the 
F.E. model. When this situation exists, either the mode shape data must be 
Interpolarted (ref. 7) or the F.E. model reduced (ref. 16). In this paper It 
will be assumed that one of these procedures has already been applied, thus 
rendering the number of degrees of freedom equal to the number of experimental 
measurement points where mode shape data Is taken. It will also be assumed 
that the experimental mode shapes are measured at the F.E. node locations. 


Sample Problem One 

Sample problem one consists of a planar cantilever beam. Two finite ele- 
ment models were used In the analysis. This first model, referred to as the 
analytical model. Is used for computing the frequencies and mode shapes that 
would normally be generated by an analytical analysis. The second model, 
referred to as the "experimental" model. Is used for simulating frequencies and 
mode shapes that would be obtained by performing an actual experimental modal 
analysis on a real beam. It Is advantageous to use simulated data In place of 
real data because the property matrices corresponding to the simulated data are 
known, whereas the property matrices for any real structure are unknown. Since 
the mass, damping, and stiffness matrix are known for the simulated experi- 
mental data, the exact error matrices can be compared to the error matrices 
generated by the equations derived In this study and the procedures can be 
evaluated. 

The analytical model Is made up of nine, equally spaced node points and 
eight connecting beam elements (fig. 1). All of the degrees of freedom are 
constrained at node 1 and every degree of freedom except for the z-dlsplacement 
and y-rotatlon are constrained at the other node points. This leaves sixteen 
active degrees of freedom for the structure. The section properties for the 
beam elements are 2.6xl0~3 for the moment of Inertia, 10x10^ for Young's 
modulus, and 2.6x10"^ for the mass density per unit length. 

The complex eigenvalue extraction solution sequence (Solution 28) of the 
MSC/NASTRAN finite element program was used to compute the free vibration 
frequencies and mode shapes for the beam. The Hess method (ref. 14) was 
selected for extracting the modes since this method Is more efficient when all 
of the modes are desired. All of the modes were Initially required for a 
complete verification of the difference matrix routines. 
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The simulated experimental model was made to differ from the F.E. model by 
adding a concentrated mass, damper, and spring to the beam. The location of 
these elements Is shown In figure 1. The values used for the elements are 
listed In table I as AH, AC, and AK. The mass, damping, and stiffness from 
the F.E. model at the same nodes and directions are also listed to give an 
Indication of the relative magnitude of the differences. NASTRAN was again 
used for computing the complex frequencies and mode shapes of the experimental 
model . 

Table II shows the comparison between the computed eigenvalues for the 
analytical model and the simulated experimental model for each of the four 
cases. All 16 of the modes were computed by NASTRAN. From table II, a com- 
parison can be made between the complex valued eigenvalues. As expected, the 
real part of the analytical eigenvalues are zero since there Is no damping 
present In the F.E. model, and the real part of the experimental eigenvalues 
are non-zero since damping is present. In general, the addition of the tip 
mass and the damper tends to lower the frequencies while the spring raises the 
frequencies. The modal damping is totally dependent on the concentrated 
damper. 

The Imaginary (frequency) part of the F.E. and experimental eigenvalues 
are plotted In figure 3 for case 1. If the eigenvalues matched exactly they 
would plot directly on the straight, 45® line. There is a small deviation from 
the straight line, but not enough to indicate any significant differences 
between the analytical and experimental models. Even If there were large 
deviations between the analytical and experimental eigenvalues, there would 
not be any way to use the results In figure 3 or table II to relate the devia- 
tions to differences in physical mass, damping, or stiffness coefficients. . 

The real components of the first four mode shapes for the analytical model 
(case 1) are plotted in figure 4. Only the translational degrees of freedom 
are plotted. Even though the first frequency has the largest deviation 
(figure 3) the first mode shape matches up very closely. The opposite occurs 
for the second and third modes where the mode shapes deviate from each other 
while the frequencies are very similar. As with the frequency plots, there Is 
no way to relate the deviation from perfect correlation in mode shape plots to 
physical differences In mass, damping, or stiffness. 

The Difference Matrix program was used to relate the differences between 
the experimental and analytical models In terms of differences In mass, damp- 
ing, and stiffness. The computer program was verified using all four cases 
and various numbers of modes as Input data. When all of the modes are Included 
as Input the only approximation In the procedure Is from the Taylor series 
truncation. As previously discussed. In a real situation all the modes would 
not be available from tests. Plots of mass, damping, and stiffness difference 
matrices for case 1, using all 16 modes are shown In figure 4. The differences 
are plotted on a grid where each intersection of a grid line corresponds to a 
location In the matrix being plotted. For example, the mass difference shown 
In the figure corresponds to the (15, 15) location In the structure's mass 
matrix. In the figure the physical differences between the analytical and 
experimental models are clearly defined. The mass difference matrix indicates 
a mass difference at degree of freedom 15 which corresponds to the transla- 
tional direction at the beam tip where the concentrated mass was added. The 
damping and stiffness errors at degree of freedom seven and one respectively, 
correspond to the locations of the concentrated damper and spring. There 
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were no other differences between the analytical and experimental models which 
Is Indicated by the flat areas In all three of the difference matrix plots. 

Even though the location of the difference Is exact, there Is some amount of 
disagreement between the actual magnitudes of the mass and stiffness, and the 
magnitudes computed by the error matrix program. The program computed a mass 
difference of 0.00011 while the actual concentrated mass was 0.0002. The 
spring magnitude was computed to be 3294, while the real spring was 5000. The 
magnitudes of the actual damper and the magnitude computed by the program were 
both 1.0. It Is not surprising that the computed mass differences was so far 
off since the mass added to the experimental model was almost as great as the 
original mass In the analytical model. Since all of the modes were Included 
In these calculations, the differences between the real values and the computed 
ones can be attributed to the higher order terms that are missing In the Taylor 
expansion. A procedure for Improving the accuracy of the magnitudes will be 
discussed later In this section. 

When less than all 16 modes are Included In the calculations the results 
deteriorate. In figure 5 results are shown for the case where only one mode 
was Included as Input Into the Difference Matrix program. The mass and damping 
difference plots do not show anything but some distributed noise. The stiff- 
ness difference plot Indicates a difference at the spring location, but the 
difference Is of the wrong magnitude. (After examining the data, the sign of 
the difference was also found to be Incorrect) Figure 6 shows difference plots 
where 10 modes are Included as Input. In this case the noise has virtually 
disappeared and the correct locations of the differences have shown up. 

A compilation of results for all four test cases are shown In figures 7 
to 9. In these figures the ratio of the computed to actual difference at the 
mass, damper, and spring location are plotted as a function of the number of 
modes used as Input data Into the Difference program. All four cases were run 
using 16, 8, 5, 3, 2, and 1 modes as Input Into the Difference program. From 
the mass difference plot (fig. 7), It Is seen that when the mass difference Is 
large (case 1), the computed difference Is only about half of the correct dif- 
ference. When the mass difference was reduced (cases 2, 3 and 4) the computed 
difference was much closer to the correct difference. If the mass difference 
Is as great or greater than the analytical mass, the location of the difference 
will be correct but the magnitude will not. From cases 2, to 4 It Is also seen 
that the computed mass difference does not change with the level of the damping 
or stiffness differences. 

Figure 8 shows the computed damping differences for the four cases. This 
plot shows that the differences are Independent of the level of damping as well 
as Independent of the magnitude of the mass and stiffness difference. Even 
when large amounts of damping are present In the structure the damping calcula- 
tions are accurate (the damping level In case 2 was close to critical). It Is 
encouraging to note that the accuracy Is Independent of the damping difference 
level, because In analytical modeling It Is the damping values that are the 
most difficult to predict. Thus for a typical structure the difference proce- 
dure would work fairly well, since the mass and stiffness differences would 
ordinarily be small, and the magnitude of the damping difference would not 
matter. 
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Figure 9 shows the effects of the various difference ratios on the com- 
puted stiffness differences. Similar to the mass calculations, the accuracy 
of the difference is dependent upon its relative magnitude. When the stiff- 
ness difference is relatively large, the computed difference is inaccurate; 
when the difference is small, the computed value is much closer to the actual 
value. Again, the computed difference is independent of the level of the dif- 
ferences in the other parameters. 

From any of the figures presented thus far it is apparent that when only 
a few modes are included the results are meaningless. When less than eight 
modes are included the results are poor, and past eight modes the results are 
good and do not improve by including more than the first eight modes. To 
determine how the number of degrees of freedom used in the model effects the 
number of modes required for good results, a new model of the cantilever beam 
was constructed using 32 degree of freedom instead of 16. 

The difference plots for this model were computed using the differences 
from case 3. The results are shown in figure 10. The difference matrices 
using 16 modes as input are shown in figure 11. From these results it is seen 
that while only 8 modes produced good results in the 16 degree of freedom 
model, at least 16 modes are needed in the 32 degree of freedom model. 

Previously, it has been shown that the accuracy of the computed differ- 
ences are dependent on the magnitude of the differences and the number of modes 
included in the calculations. In an attempt to improve the accuracy an itera- 
tive procedure was implemented (fig. 12). In this procedure the differences 
computed by the Difference program are accumulated from all previous Iterations 
and are then added to the mass, damping, and stiffness matrices for the analyt- 
ical model. The updated analytical model is then used to compute a new set of 
differences for the next iteration. 

The iterative procedure was tested using the differences from test case 3 
and the sixteen degree of freedom model. The results for three iterations are 
shown in figure 13. Without iterating, it was shown that when all sixteen 
modes are included in the calculations the results are very good. After iter- 
ating only twice, the results are exact. The same is also true when only eight 
modes are used. In general, when less than eight modes are used, the accuracy 
of the computed differences are not improved significantly by iterating. When 
only a few modes are included, the accuracy is not improved at all. The 
advantage of using the iteration procedure is that when an adequate number of 
modes are used the results will converge to the exact values regardless of the 
magnitude of the differences. The limitation of the iterative process is that 
it does not reduce the number of modes required for good results. 


SAMPLE PROBLEM TWO 

The second sample problem consists of a planar, simply supported beam. 
The finite element model of this problem is made up of nine node points and 
eight connecting beam elements (fig. 14). All of the degrees of freedom are 
constrained, except for the y-rotatlons at nodes one to nine, and the 
z-translations at nodes two through eight. There are sixteen degrees of 
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freedom for this problem. The same section properties that were used of for 
the first sample problem are also used here. The difference matrix plots for 
this problem were generated using the Iteration scheme shown In figure 16. 

The "experimental" model was made to differ from the analytical model by 
adding three concentrated springs and seven dampers to the beam model. The 
locations and properties for these elements are shown In the figure. This 
sample problem differs from the first one In that the differences In the first 
problem were limited to a single mass, damper, and spring, while In this 
problem there are several springs and dampers at every node. Also, the level 
of damping Is much less In this problem than In the previous one. 

A comparison of the eigenvalues for the second sample problem Is shown In 
table III. From this comparison It Is seen that the major differences between 
the analytical and experimental models are In the first frequency and the modal 
damping In the first seven modes. Beyond the seventh mode there are not any 
differences between the analytical and experimental eigenvalues. The first 
frequency is higher for the experimental model because of the additional stiff- 
ness from the three springs. The modal damping Is different because the exper- 
imental model has the seven translational dampers while the analytical model 
does not have any damping. It Is understandable that there is no modal damping 
In the higher modes for the experimental model since the higher modes are 
dominated by rotations and the dampers only act in the translation direction. 

The computed damping and stiffness difference matrices using the first 
mode only as Input Into the Difference program are shown In figure 15. From 
these plots It Is Impossible to Identify any of the differences that actually 
exist between the analytical and experimental models. When only one mode was 
used for the first sample problem calculations the differences could not be 
Identified either. When the number of modes was increased to four (fig. 16) 
the correct differences were reasonably apparent In the difference plots. For 
the stiffness matrix plot, differences appear at degree of freedom 4, 6, and 8 
which corresponds to the locations of the three springs that were added to the 
experimental model. In addition to the differences at these degrees of free- 
dom, differences also appear at some of the other degrees of freedom. These 
differences do not actually exist In the models and would not appear If more 
modes were used as Input. Many of these "extra" differences can be eliminated 
by examining the possible coupling that mayexist In the analytical mode. For 
example, node two and six are not connected to each other so degrees of freedom 
2 (z-translatlon, node two) and degree of freedom 10 ( z-translatlon, node six) 
are uncoupled which allows for location (2,10) and (10,2) In the difference 
matrix to be set to zero. The same logic can be used to eliminate some of the 
other unobtainable differences appearing in the plots. 

In figure 17, the difference plots using four modes are recreated, except 
that the differences at uncoupled degrees of freedom are set to zero. In these 
plots the correct differences are even more evident although some differences 
continue to appear where there are not any true differences. There does not 
appear to be any way to eliminate these "extra" differences except by using 
more modes In the calculations. When the number of modes Is Increased to six 
(fig. 18), both the computed stiffness and damping difference matrices are very 
accurate. In addition, no significant differences appear where they do not 
actually exist In the models. In figure 19, plots are shown where all sixteen 
modes are Included In the computations, and as expected, the results are 
almost exact. 


11 



CONCLUSION 


A general procedure for identifying and quantifying the differences 
between F.E. models and experimental data has been developed and demonstrated 
with simulated experimental data. The differences, which can be computed for 
linear, viscously damped structures, are presented In terms of mass, damping, 
and stiffness coefficients. Since the differences are computed In terms of 
mass, damping, and stiffness coefficients, possible modeling problems can be 
Identified In the F.E. or analytical model. 

From data generated for a damped cantilever beam and a damped simply sup- 
ported beam. It was determined that the accuracy of the computed differences 
Increases as the number of experimentally measured modes Included In the cal- 
culations Is Increased. When the number of experimental modes is at least 
equal to the number of translational degrees of freedom both the location and 
magnitude of the differences can be computed very accurately. When the number 
of modes Is less than this amount the location of the differences may be 
determined even though their magnitudes will be under estimated. When too few 
modes are available neither the location or the magnitudes of the differences 
can be identified. 

In practice. It will be required to measure the experimental frequencies 
and mode shapes very accurately before the differences can be attributed to 
shortcomings In the analytical model. If the experimental data Is not precise, 
the computed differences can still provide considerable Insight Into the pos- 
sible locations of deficlences. The difference Is that the deficiencies may be 
In the experiment and some judgement will be required to decide whether to 
modify the experiment or the analytical model. 
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TABLE I. - TEST CASES FOR 
SAMPLE PROBLEM ONE 


Case 

AH/H^ 

AC/Cb 

Ak/kC 


2x10“* 

1.0 

5000 


2.6x10“* 

0 

78000 


5x10“^ 

1.8 

2500 


2.6x1 0“* 

0 

78000 


5x10“® 

1.8 

2500 

3 

2.6x10“* 

0 

78000 

4 

5x10“® 

1 .0 

25000 

2.6x10“* 

0 

78000 


*Rat1o of mass A to F.E. 

mass at node 9. 

^Rat1o of damping A to F.E. 

damping at node 5. 

CRatIo of stiffness A to 
F.E. stiffness at node 2. 


TABLE II. - COMPARISON OF ANALYTICAL AND "EXPERIMENTAL" EIGENVALUES FOR SAMPLE PROBLEM ONE 


Mode 

Analytical 

eigenvalue 

Case 1 

Case 2 

Case 3 

Case 4 

Experimental 

eigenvalue 

Experimental 

eigenvalue 

Experimental 

eigenvalue 

Experimental 

eigenvalue 

1 

0 


1371 

-44.3 


1301 

-123 

♦ 

nil 

-53 

♦ 

1331 

-45 

♦ 

1501 

2 

0 

♦ 

84D1 

-247 


7611 

-438 


5811 

-245 


7821 

• 258 


8741 

3 

0 

* 

2 3121 

-179 


2 2711 

-12 


2 3141 

-74 


2 3171 

-31 


2 5251 

4 

0 


4 4471 

-203.9 


4 3811 

-378 


4 3591 

-21.4 


4 4251 

-183 


4 7941 

S 

0 


7 1971 

-28.2 


7 1781 

-30 


7 2031 

-17 


7 2071 

-59 


7 6921 

6 

0 

* 

10 4271 

-195 


10 4281 

-360 


10 3871 

-202 


10 4301 

-158 


10 9721 

7 

0 


13 7571 

-247 


13 7941 

-35 


13 7761 

-19.6 


13 7801 

-55 


14 2351 

6 

0 


16 3711 

-197 


16 3691 

-359 


16 3041 

-199 


16 3591 

-168 


16 5721 

9 

0 

+ 

63 0171 

0 

♦ 

62 9221 

0 


62 9821 

0 


62 9821 

0 


62 9821 

10 

0 


75 7271 

0 


75 7111 

0 


75 7211 

0 


75 7211 

0 


75 7211 

11 

0 


84 0331 

-3 


83 9441 

-5 


840191 

-3 


84 0191 

-3 


84 0191 

12 

0 


94 0261 

0 


93 9741 

0 


94 0071 

0 


94 0071 

0 


94 0071 

13 

0 


103 9281 

-2 


103 8751 

-4 


103 9091 

-2 


103 9091 

-2 


103 9091 

14 

0 

* 

112 5991 

0 


112 5591 

0 


112 5841 

0 


112 5841 

0 


112 5841 

1$ 

0 


119 2781 

0 


119 2561 

-1 


119 2701 

0 


119 2701 

0 


119 2701 

16 

0 


123 4711 

0 


123 4661 

0 


123 4701 

0 


123 4691 

0 


123 4701 


TABLE III. - COMPARISON OF 
ANALYTICAL AND 
■EXPERIMENTAL" 
EIGENVALUES 
FOR SAMPLE 
PROBLEM 


TWO 


Mode 

Analytical 

eigenvalue 

Experimental 

eigenvalue 

1 

0 t 3861 

-48 + 4881 

2 

0 + 1 5411 

-48 + 1 5601 

3 

0+3 4551 

-48 + 3 4661 

4 

0+6 0931 

-48 + 6 0951 

5 

0+9 3491 

-47 + 9 3521 

6 

0+12 9171 

-47 + 12 9191 

7 

0 + 16 0311 

-48 + 16 0331 

8 

0+62 6681 

0+62 6681 

9 

0+62 9301 

0+62 9301 

10 

0 + 76 6171 

0 + 76 6171 

11 

0+85 5771 

0+85 5771 

12 

0+95 5621 

0+95 5621 

13 

0 + 105 1191 

0 + 105 1191 

14 

0 + 113 3471 

0 + 113 3471 

15 

0 + 119 6331 

0 + 119 6331 

16 

0 + 123 5631 

0 + 123 5631 



















Relative displacement 



Figure 3. - Comparison of mode shapes for sample problem one, case 1. 





ABS difference 



A8S differenc* 



(al Mass dlllwence matrix. 



(M Stillness dillerence matrix. (c) Damping difference matrix. 

Figure 6. - Computed differences for sample problem one, case 1, using 10 modes. 


Computed difference/exact difference Computed difference/exact difference 



Figure 7. - Mass difference at concentrated mass location (sample problem one). (*See 
table I for case definition.) 



Figure 8. - Damping difference at concentrated damper location (sample problem one). 
("See table I for case definition. ) 
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Number of modes 


(c) Stiffness difference. 

Figure 10. - Computed differences for 32 D.O.F. model. 


ABS difference 



(b) Damping difference matrix. (c) Mass difference matrix. 

Figure U. • Computed differences for 32 O.O.F. model, using 16 modes. 


Computed difference/exact difference 



Figure 12. - iteration scheme. 
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(a) Mass difference. (b) Damping difference. 


(c) Stiffness difference. 


Figure 13. - improvement in computed differences with iterations (sampie problem one, 
case 3). 
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(a) Analytical model. 


I r f 



0.05 

(Typical) 


(b) "Experimental" model. 
Figure 14. - Sample problem two. 




(a) Stiffness difference matrix. (b) Damping difference matrix. 

Figure 15. - Computed differences for sample problem two, using 1 mode. 
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(a) 


Stiffness difference matrix. 

Figure 16. - Computed differences for 


(b) Damping difference matrix, 
sample probiem two, using 4 modes. 
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ABS difference 







(a) Stiffness difference matrix. (b) Damping difference matrix. 

Figure 17. - Computed differences for sample problem two, using 4 modes and coupling constraints. 





(a) Stiffness difference matrix. 


(b) Damping difference matrix. 


Figure 18. - Computed differences for sample problem two, using 6 modes. 




ABS difference 



(a) Stiffness difference matrix. (b) Damping difference matrix. 

Figure 19. - Computed differences for sample problem two, using 16 modes. 
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